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The successive perturbative estimates of the pressure of QCD at high temperature T show 
no sign of convergence, unless the coupling constant g is unrealistically small. Exploiting 
known results of an effective field theory which separates hard (order 2ttT) and soft (order 
gT) contributions, we explore the accuracy of simple resummations which at a given loop 
order systematically treat hard contributions strictly perturbatively, but soft contributions 
without truncations. This turns out to improve significantly the two-loop and the three-loop 
results in that both remain below the ideal-gas value, and the degree of renormalization scale 
dependence decreases as one goes from two to three loop order, whereas it increases in the 
conventional perturbative results. Including the four-loop logarithms recently obtained by 
Kajantie et al., we find that this trend continues and that with a particular sublogarithmic 
constant the untruncated four- loop result is close to the three-loop result, which itself agrees 
well with available lattice results down to temperatures of about 2.5T C . We also investigate 
the possibility of optimization by using a variational ( "screened" ) perturbation theory in the 
effective theory. At two loops, this gives a result below the ideal gas value, and also closer to 
lattice results than the recent two-loop hard-thermal-loop-screened result of Andersen et al. 
While at three-loop order the gap equation of dimensionally reduced screened perturbation 
theory does not have a solution in QCD, this is remedied upon inclusion of the four-loop 
logarithms. 



I. INTRODUCTION 

One could expect weak coupling calculations to lead to reasonable estimates of the QCD free 
energy at high temperature T, a regime where indeed the gauge coupling becomes small because 
of asymptotic freedom. But explicit perturbative calculations, which have been pushed in recent 
years up to the order g 5 jH S S IE IE 0] do not exhibit any sign of convergence, as depicted 
in Fig. ^ they rather show increasing ambiguities due to the dependence on the renormalization 
point, signalling a complete loss of predictive power. 

Various mathematical extrapolation techniques have been tried, such as Pade approximants 
|E 0] and Borel resummation • The resulting expressions are indeed smooth functions of the 
coupling, better behaved than polynomial approximations truncated at order <? 5 or lower, with a 
weaker dependence on the renormalization scale. However, while these methods do improve the 
situation somewhat, it is fair to say that they offer little physical insight on the source of the 
difficulty. 

Recognizing that an important effect of thermal fluctuations is to give a mass to the excitations, 
thereby screening long range interactions, it has been suggested to incorporate such screening effects 
in the tree-level Lagrangian, and correct for double counting by adding suitable counterterms to the 
interactions. Such a scheme has been implemented with some success in scalar field theory under 
the name of screened perturbation theory (SPT) 0, ^ l3.ll3|- It has been extended to QC D by 



including at tree level the entire non-bear Ls^TJ^ ^ard thermal loops Q Q, Q 0, 
which is referred to as HTL perturbation theory (HTLPT). 

A different approach, motivated physically by the success of the quasiparticle picture, is based 
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FIG. 1: Strictly perturbative results for the thermal pressure of pure glue QCD normalized to the ideal-gas 
value, as a function of T/T c (assuming T c /Ky^ = 1.14). The various gray bands bounded by differently 
dashed lines show the perturbative results to order g 2 , g 3 , g 4 , and g 5 , with MS renormalization point p, 
varied between ttT and AirT . The thick dark-grey line shows the continuum-extrapolated lattice results from 
Ref. [2^; the lighter one behind that of a lattice calculation using an RG-improved action |27j . 



on the so-called ^-derivable approximations lid, |20|, UM ( see also Peshier j22|). This approach 
takes advantage of remarkable simplifications which occur in the calculation of the entropy at two 
loop in the skeleton expansion. Together with further approximations 1 for the self-energies based 
on hard thermal loops, this led to results for the thermodynamical functions which are consistent 
with lattice calculations for temperatures above 2.5 T c . 

This paper reconsiders the known results up to order g 5 in the light of the simple observation 
that the accuracy of perturbation theory is not the same at all momentum scales. Thus, while 
perturbation theory at the scale T is an expansion in powers of g 2 , the perturbation theory at 
the scale gT is an expansion in powers of g, and is therefore less accurate. It is when they are 
treated strictly perturbatively that the soft contributions turn out to completely spoil apparent 
convergence, as exemplified by the (soft) contribution of order g 3 which leads to a pressure exceeding 
the ideal gas value. The main idea that we want to pursue is to decouple approximations in the hard 
and the soft sectors: in the hard sector, we shall use perturbation theory since it is accurate; in the 
soft sector we shall go (minimally) beyond perturbation theory, with the remarkable result that the 
apparent convergence of perturbative QCD at high temperatures is dramatically improved. This, 
in our opinion, lends support to the more ambitious attempts to reorganize perturbation theory 
by novel resummation techniques. 

The tools to deal with various momentum scales in field theory are well developed and involve 
the construction of effective theories. In fact, the perturbative results through order o 5 , first 
calculated in Refs. 0, E|, have been rederived and confirmed by Braaten and Nieto Q by an 
elegant and efficient effective-field-theory method which separates contributions from hard, soft, 
and supersoft momentum scales, 27rT, gT, and g 2 T, respectively. In Euclidean space, the only soft 
modes are static ones, and the effective field theories for these are therefore dimensionally reduced 



1 A full "^-derivable approximation has been worked out successfully in scalar field theory |2al . for which recently 
also the question of renormalizability could be answered affirmatively in Refs. [23. 12^1 . 
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|2ft I29L |3CJ to three-dimensional ones. The dimensionally reduced theory consists of massive 
adjoint scalar fields Aq and massless three-dimensional Yang-Mills fields (with coupling gg) and, 
following Braaten and Nieto, we shall refer to it as electrostatic QCD (EQCD) in the following. 
The corresponding effective three-dimensional Lagrangian is 

Ceqcd = \f«F* + ^(DiA ) a (DiA ) a + \m\A%A% + \\ E {A%A%f + 5C E qcd (1) 

where the parameters are determined perturbatively by matching 7j. In lowest order we have: 

m 2 E = (1 + N f /6)g 2 T 2 , g\ = g 2 T, (2) 

and 0] 

In fact, A^; starts to contribute to the pressure only at order g 6 . For this reason we shall ignore this 
particular coupling in most of the following, as well as all the other vertices contained in 5 Ceqcd- 

The thermal pressure of the 4-dimensional theory can be decomposed into contributions from the 
hard modes ~ T, calculable by standard perturbation theory, and soft contributions governed by Q 
which involves both perturbatively calculable contributions up to order g 5 T and nonperturbative 
ones coming from the fact that the effective theory for the modes Ai{x) is a confining theory. 
However the latter magnetic contributions start at order g e , and cannot a priori be made responsible 
for the poor apparent convergence that is seen up to order g 5 . 

We shall then focus in this paper on the soft contributions, i.e. those coming from the momentum 
scale gT. In principle these can be calculated using perturbation theory, but as clear from (^Q) and 
(J2J the corresponding expansion parameter is g E /rriE ~ g, so that the perturbative expansion 
converges only slowly, more slowly than the perturbative expansion in the hard sector. As a 
minimal step towards a nonperturbative treatment of the soft sector, we shall perform a simple 
loop expansion of keeping the parameters tue and gE as given in terms of g by the matching 
conditions, and not expanding them out in powers of g in the final result. As we shall see, this 
simple method leads to a significant improvement over the strict perturbative results. We also 
consider the effect of including the four-loop logarithms recently obtained by Kajantie et al., and 
again find that strict perturbation theory has large scale dependences which are drastically reduced 
when keeping soft contributions untruncated. For a particular choice of the sublogarithmic constant 
the untruncated four- loop result is moreover close to the three- loop result, which itself agrees well 
with available lattice results down to temperatures of about 2.5T C . 

In Sect. Mil we consider a simple variational improvement of perturbation theory in the form 
of dimensionally reduced screened perturbation theory (DRSPT). This turns out to be much sim- 
pler than HTLPT, while also allowing to resum screening masses in a gauge invariant manner. 
At two-loop order it leads to a result significantly closer to lattice data than two-loop HTLPT, 



which suggests that the partial failure of the latter as observed in [17|, [18| is due to spurious hard 
contributions. At three-loop order, the gap equation of DRSPT has no real solution, but this is 
remedied upon inclusion of the four-loop logarithms. 

In the final section we summarize and discuss our results and try to put them into perspective 
with other techniques. 



II. DIMENSIONAL REDUCTION BEYOND STRICT PERTURBATION THEORY 

In the following we adhere to Ref. [3] in the treatment of the dimensionally reduced effective 
theory responsible for the contributions from the scale gT, but we deviate in that we shall not 
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treat the soft sector strictly perturbatively. We shall organize our presentation by considering 
the successive approximations obtained by expanding the contributions of the hard modes to the 
pressure in powers of a s = g 2 /(4-7r). Each order in this expansion defines also the accuracy with 
which the parameters of the effective Lagrangian are determined through perturbative matching 
conditions. However, at each level of approximation we shall consider the effect of treating the 
contributions of the soft sector more completely by refraining from the specific truncations usually 
employed in a strictly perturbative expansion. 



A. One-loop order 

In massless QCD at one-loop order, the only contribution to the thermal pressure (identical to 
minus the free energy) is coming predominantly from hard momenta. Introducing a momentum 
cutoff A^; to separate the hard scale 2irT from the scale gT, one gets a contribution from the 
soft sector proportional to TA E , compensating a similar term in the interaction- free one-loop 
contribution from hard modes. One ends up with the standard ideal-gas result Pq = 7r 2 T 4 (8/45 + 
7Nf/60) for Nf quarks and 3 colors. 

Ref. Q avoids the introduction of momentum cutoffs by using dimensional regularization for the 
purpose of both ultraviolet and infrared regularizations. Doing so, the one-loop result exclusively 
comes from the hard modes: At this level of approximation where all interactions are neglected, 
the soft (zero) modes are to be taken as massless, and their contribution vanishes in dimensional 
regularization where scaleless integrals are set to zero. 



B. Two-loop order 

The two- loop, i.e. order a s , contribution of hard modes to the pressure is still independent of 
a cutoff Ae, if this is handled by dimensional regularization. It reads [lj 

^a! d = -f(l + ^/H^. (4) 

At this level of approximation, the soft modes described by EQCD are massive, the mass 
being given by the leading order Debye mass The one- loop contribution to the pressure from 
EQCD is now, on dimensional grounds, proportional to m? E times an overall factor of T. If only 
this contribution is added on to (@J, the result is the ill-behaved strictly-perturbative result to 
order g 3 displayed in Fig. ^ However, in the soft sector, we need not restrict ourselves to this 
one-loop approximation, but rather treat more completely the interactions which are present in 
C-eqcd- Thus, with the parameters of EQCD determined by matching with the perturbative 
calculation at the present level of accuracy, we shall consider also the two-loop contributions from 
the dimensionally reduced effective theory. This yields the following contribution to the pressure 

p(l)+(2) _ T ,(i)+(2) 

EQCD _ 1 JM \°) 

with f$ +{2) given by 

f (i)+(2) _ 2 3 3 (I Ae ^ „\ 22_ 



ft™ = -3-4 + W ( 7 + 41n ^ + 3 ) 9E m E + Sf E , (6) 



where A^ is the scale of dimensional regularization in the soft sector, which may be loosely associ- 
ated with the separation scale between the soft and the hard momenta. Choosing the counterterm 
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SJe by minimal subtraction, 

5 fE = -g^29E m E-: ( 7 ) 

leaves behind a dependence on A^. Eventually, ln(A£;) has to combine with a matching logarithm 
arising from the hard scales for which Ae provides the infrared cutoff. This in fact happens after 
a careful perturbative matching at three-loop order (see eq. ©)• 

Together @ and (JSJ) are accurate through order g 4 log(l/g), with an error of order g 4 as to be 
expected from a two-loop calculation. The coefficient of the g 4 log(l/g) term is correct provided 2 
Ae is set to a constant times T. From its role in dimensional reduction, it should be smaller than 
2irT but larger than gT. However, as remarked in Ref. Q, the introduction of cutoffs through 
dimensional regularization leaves their relationship to momentum cutoffs undetermined, and there 
could be a different relationship between the scale of dimensional regularization and effective mo- 
mentum cutoffs depending on whether the latter act as IR or as UV cutoffs. Ref. even found 
that the scale Ae should be chosen larger than the UV scale fl ~ 2irT in order to have opti- 
mal convergence of the hard contributions. For simplicity, and to facilitate the comparison with 
other calculations where a similar identification is made, we shall put A^; = fl in the following, 
i.e. introduce only one scale for dimensional regularization, and refrain from modifying it by hand 
depending on whether the various logarithms can be identified as arising from regularization in the 
IR or in the UV. 3 
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FIG. 2: Two-loop pressure in pure-glue QCD with untruncated EQCD contributions when fl is varied 
between ttT and 4ttT (broad gray band) in comparison with the lattice result from Ref. (thick dark- 
grey curve). The narrow darker-grey band above the former is the result of 2-loop DRSPT considered in 
Sect. IIII Al its lower boundary corresponds to the extremal value when varying ft. 



2 If Ae is chosen to be parametrically smaller than 2tyT by multiplying it by a fractional power g c with < c < 1, 
then the coefficient of g 4 ln(l/g) is wrong by a factor of c. 

3 The ambiguity of the choice of Ae can alternatively be understood as arising from the freedom to renormalize 
the effective 3-dimensional theory differently than by minimal subtraction. Because in the next subsection we 
shall compare with HTLPT [Til ll(l ITtI Hjj| where only minimal subtraction of additional divergences has been 
considered, we stick to minimal subtraction in the following. It should be kept in mind, however, that there is a 
source of additional renormalization scheme dependence. 
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In Fig. [21 we give the numerical evaluation of the full 2-loop result obtained as indicated above, 
for pure-glue QCD and T between T c and 5T C in analogy to the strictly perturbative results of 
Fig. ^ We always use the standard 4 perturbative solution to the two-loop renormalization group 
equation for a s assuming T c /Aj^g = 1.14. The UV renormalization scale ft is varied about a central 
value 2ttT by a factor of 2, and it should be kept in mind that this variation also traces some of 
the ambiguity in choosing Ar. The resulting total scale dependence is comparable to the scale 
dependence of the perturbative order-g 4 result. But in contrast to both the result to order g 3 , 
to which it is perturbatively equivalent, and the order-g result, the untruncated 2-loop result of 
dimensional reduction remains below the ideal-gas, and thus has overlap with the lattice results, 
which the former do not have. 

Remarkably, the partial inclusion of order-g 4 effects arising from a two-loop evaluation of EQCD 
is superior to a complete order-g 4 evaluation in strict perturbation theory. Although the former 
has an uncancelled A^ dependence in addition to the normal /2-dependence, numerically the total 
scale dependence is not worse but even slightly better than that of the perturbative order-g 4 result. 

If we had not put A^ to zero on grounds that it starts contributing only at order g e , we would 
have obtained the additional 2-loop term 



f (2) 
J M 



5 

^m 2 E X E . (8) 



Inserting the leading-order value of As, Eq. ©, this contribution is not only of order g 6 , but it 
is also numerically quite small in comparison with the other two-loop contributions © even when 
a ~ 1. 



C. Three-loop order 

The three- loop contribution of the hard modes to the pressure is no longer Ag-independent, 
because it has to be matched with the minimally subtracted, and thus A^-dependent, EQCD 
theory. This has been calculated in [?] as 5 



fid = ^ T4 { 244 - 9 + 17 - 2m f ~ 0A15N f + 135 f 1 + ln 



a 



2 



1 + ^K 1 -#^J ta 5SrJlTJ- (9) 

Similarly, the mass parameter of EQCD can be obtained by a matching calculation of two-loop 
self-energies as 0] 



5 + 22 7 + 22 In ; ' 



(i_ 8 i n2 + 7 7 + 71n-4^) fl-2 7 -21n4;ll) 



4vrT 

(10) 



As for the three-loop contribution from the soft sector, this is given by the finite and thus 
Ag-independent expression calculated in Ref. 7] (neglecting A_E-contributions now) 

( 3) 9/89 11 1 2 \ 4 

tW = n m2H 7T OpTTT-p. 11 

JM 8?r 3 V24 6 6 J E E y 1 



4 This is practically indistinguishable from the full 2-loop solution as soon as T > 2T C (see Appendix of Ref. H^])- 

5 In Eq. (54) of Ref. the terms corresponding to the second and fourth term on the r.h.s. of Eq. Q have a 
different sign due to a typographical error. 
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In a strictly perturbative evaluation which drops all terms of order g°, the sum of the hard and 
soft contributions are Ag-independent as they should be. That is, at order g 4 , the term oc lnA^; 
in Q cancels against the corresponding one in ©. These perturbative contributions, evaluated 
numerically, give the result marked 

V in Fig. [U or Fig. EH 

However, our aim is to go beyond such perturbative results, and as a simple approximation 
in this direction, we consider keeping the soft contributions @ and Ijlljl untruncated when the 
perturbatively determined m E is inserted. This then corresponds to a selective summation of 
higher-order effects that may help improve the convergence of perturbation theory, although these 
higher order terms are Ag-dependent. 6 
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FIG. 3: Three-loop pressure in pure-glue QCD with untruncated EQCD contributions when g, is varied 
between irT and AttT (medium-gray band); the dotted lines indicate the position of this band when only 
the leading-order result for toe is used. The broad light-gray band underneath is the strictly perturbative 
result to order g 5 with the same scale variations. The full line gives the result upon extremalization (PMS) 
with respect to p. (which does not have solutions below ~ 1.3T C ); the dash-dotted line corresponds to fastest 
apparent convergence (FAC) in m%, which sets g rs 1.79ttT. 



Setting Ae = fl the ambiguity in choosing A^ contributes to the scale dependence and is 
included in our error bands as g is varied around its central value. 7 The result is shown in Fig. |3J 
Despite the incomplete cancellation of Ae in the untruncated evaluation of the three- loop result, 
we observe a considerable reduction of the total scale dependence compared to strict perturbation 
theory. 

Furthermore, compared to the full two-loop result, the three-loop result stays within the (rather 
large) uncertainties of the former. While the uncertainties remain sizeable even at three-loop order, 
the overall picture is a remarkable improvement over strict perturbation theory, where the results 
jump about the ideal-gas value and the renormalization scale dependence increases steadily with 
the highest power of g reached. 



6 Specifically, they lead to a g 6 ln(g) contribution with the constant under the log carrying the ambiguity in Ae if 
the latter is proportional to T. Numerically, however, this g 6 \n(g) contribution is completely negligible compared 
to the recently determined |33« perturbative g 6 ln(g) contribution appearing at 4-loop order. 

7 These error bands would of course be widened by an independent variation, so the scheme dependences displayed 
in our figures are certainly underestimated to some extent. 
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FIG. 4: PMS-extremalized full-three-loop pressure in QCD with Nf = (full line), Nf = 2 (dashed line), 
and Afy = 3 (dash-dotted line) in comparison with the estimated continuum extrapolation of QCD with 2 
light quark flavors of Ref . [34| • 

It is interesting to observe that very similar results are obtained if only the lowest order Debye 
mass is used in this calculation, in place of the full order g 4 expression (jlUj) . This is displayed in 
Fig. |31 by the dotted lines. It indicates that what matters here is the accuracy with which one treats 
the soft sector, more than the accuracy with which the parameters of Ceqcd are determined. 

Whereas in the untruncated two-loop result the scale dependence is monotonic and does not 
allow for its elimination by a principle of minimal sensitivity (PMS), such an elimination turns out 
to be possible in the three- loop result 8 . Choosing a s {fi) according to the 2-loop renormalization 
group equation, we find an extremum of the untruncated three-loop result as a function of Jx. As 
shown in Fig. |HJ the corresponding pressure values are in fact remarkably close to the lattice data 
for T > 3T C . 

The extremum is only a local one with respect to p,. For large temperatures T > 10T C it occurs 
at /2 ~ 2ttT, which following Ref. 7] we have taken as central value because it is the spacing of the 
Matsubara frequencies. For smaller temperatures, the required value of Jl increases and exceeds 
AttT below T ~ 3T C , where the lattice data start to deviate from the three-loop pressure. For still 
smaller temperatures T < 2T C the required Jl becomes unreasonably large, and finally for T < 1.3T C 
the local extremum disappears completely. (The strictly perturbative result, on the other hand, 
has a monotonic, i.e. run-away, scale dependence for all T.) 

For completeness, we also give the numerical results obtained by including Nf = 2 and 3 
massless quark flavors. The scale selected by the minimal sensitivity turns out to be slightly higher 
than at Nf = 0, it also becomes large for T < 2T C , but the extremum exists down to T c now. There 
exist no reliable continuum extrapolated lattice data for this case yet, but in Ref. 34] an estimated 
continuum extrapolation has been given for Nf = 2 light quark flavors. This is compared with 
the extremalized full 3-loop results in Fig. 0] The fact that the cases Nf = 2 and 3 are nearly 
degenerate (when normalized to the respective ideal-gas values and plotted as a function of T/T c 
with the respective critical temperature [assumed to be 1.14Aj^g]) is consistent with lattice results 
[34| , and is very similar to the pattern observed in the "NLA" results of Ref. .20] for the entropy. 



This is possible in fact only if the complete order-g 4 expression for the Debye mass is used. 
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D. Four-loop order 



In an impressive four-loop calculation, the authors of Ref. |3Jj|] have recently determined the 
last coefficient in the perturbative expansion of thermal pressure of QCD that can be computed 
analytically. At four-loop order in the effective theory Q there appear two logarithmic terms 
whose coefficents have been obtained as |3c 



p (4) /T 



N, 



(4vr) 4 



43 
12 



1577T 

768 



21 



1 AE + 

ln-n- + 
9e 



43 
T 



491vr 2 
768 



In — — + c 

m E 



(12) 



where the first logarithm is in fact from the magnetostatic sector. To obtain the complete g G Ing- 
contribution in the QCD pressure one now also needs g\ to order g 4 , given in Ref. [HH] as (for 
N f = 0) 



9% 



4vra s 1 + 



Air L 



22(ln 



4vrT 



+ 7) + l 



(13) 



The constant c in (|12|) . however, is strictly nonperturbative and can in principle be determined 
by lattice simulations of the dimensionally reduced theory, but require even more (also analytical) 
work 33]. 

Also, the four-loop contribution to Phard is not yet known, but all the terms proportional to g 6 
and involving explicit logarithms of /2 or A E are determined by the p, and independence of the 
total pressure to order g 6 , and have been given explicitly in Ref. [33j — only a constant times g 6 
thus remains undetermined. Equating again A E with ft, the four-loop contribution to be added to 
the above three-loop one when evaluated with (|13|) can be written as (again for Nt = only) 



P (4) = Z!Lt 
45 



vr 2 . f 21945 



16 



In 2 



JJ_ 
2ttT 



+ 2676. 4 In 



2vrT 



a. 



IT 



+ 



27 

327T 4 



"43 


157vr 2 " 


, T 
ln^r + 

9% 


"43 


491vr 2 " 


12 ~ 


768 


T 


768 



In— + 7.575) <il T. 

m E 



(14) 



where for easier comparison with Ref. [331 ] we have collected all undetermined constants to this 
order in a new constant 5, chosen such that the g E T contribution in ()14l) . when expanded in g, is 
proportional to (ln(l/g) +5). 

The unknown coefficient 5 of course leaves the numerical outcome completely open until the 
required analytical and numerical calculations will have been performed. However, it has been 
observed in Ref. 33] that this coefficient could well be such that the perturbative result follows 
closely the 4-d lattice results. In obtaining numerical results, Ref. [23j in fact used a particular 
optimized renormalization scheme, introduced in Eill . which also involves keeping the parameters 
of the effective theory unexpanded. (Refs. pifll l3fij ] also mentioned that this reduces the scale 
dependence.) 

In Fig. [21 we present our numerical results obtained by adding the 4- loop terms (|14|) to the full 
3- loop result of the previous section, but now evaluated with ()13j) in all terms involving g E , for 
the possibilities 5 = 0, 1/3, 2/3. This differs from [33j in that we treat hard contributions strictly 
perturbatively and only soft ones without truncations. Furthermore, we use the perturbative 2- 
loop running coupling 9 and vary Jx about 2nT by a factor of 2 as above. The results, which are 



Ref. |-'~)3l found S ~ 0.7 to give results which agree well with the 4-d lattice results. The difference to our central 
value of 1/3 is mainly due to the fact that we used 2-loop rather than 1-loop running coupling. Like Ref. [33| we 
included also the contribution JHJ, which is however fairly small. 
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FIG. 5: Four-loop pressure in pure-glue QCD including the recently determined <? 6 ln(l/g) contribution of 
[13 together with three values for the undetermined constant 6 in Eq. (|14(l when evaluated fully with p, 
varied between ttT and 4ttT (medium-gray bands). The broad light-gray band underneath is the strictly 
perturbative result to order g 6 corresponding to the central value 6 — 1/3. The full line gives the untruncated 
result with 6=1/3 extremalized with respect to /2 (which does not have solutions below ~ 1.9T C ); the dash- 
dotted line corresponds to fastest apparent convergence (FAC) in m\, which sets p, « 1.79ttT. 



full 4-loop, (5 = 0, 1/3, 2/3 




displayed by the medium-gray areas, show a remarkably small scale variation. By contrast, the 
strictly perturbative evaluation (given for 5 = 1/3 only) shows an increased scale dependence when 
compared to the 3- loop results. 

Also like in the previous 3-loop case, we find that the untruncated result has again a nonmono- 
tonic scale dependence which makes it possible to fix the scale by PMS. 10 The result is again close 
to the FAC choice considered previously as well. The strictly perturbative result on the other hand, 
has a run-away scale dependence for almost all T. 



III. SCREENED PERTURBATION THEORY IN THE SOFT SECTOR 

In massless scalar field theories, which have a poorly convergent perturbative series for the 



11 



thermal pressure similar to what is found in QCD, Karsch, Patkos and Petretcky and others 
EES have proposed to ameliorate the situation by a variationally improved perturbation theory 
which uses a simple mass term as variational parameter. In this so-called screened perturbation 
theory (SPT) the mass term is part of the tree-level Lagrangian as well as the interactions, where 
it is counted like a one-loop counterterm. 

This approach has been extended to QCD by Andersen et al. 2> 3 12 > 3] by using the hard- 
thermal-loop (HTL) action in place of a simple mass term, turning its prefactor, which is propor- 
tional to the Debye mass squared, into a variational parameter. This HTL-screened perturbation 
theory (HTLPT) has been recently carried through to two- loop order in QCD in Refs. |l7Lll«|. The 
result is perturbatively correct to order g 4 \n(g) and has been found to give rather stable results 
which are smaller than the ideal-gas pressure, but significantly above the lattice results. 



10 This has been observed also before in Ref. |3t|. but using one- loop running and the particular parametrization of 
Ref. E3. 
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In this section we investigate the possibility of improving the soft contributions considered above 
to two and three loop order by a dimensionally reduced screened perturbation theory (DRSPT) 
for EQCD defined by trivially rewriting the EQCD Lagrangian according to 

£eqcd = + i(A4>)°(A4>) a + + 6m 2 )A%A a - hm 2 A a A a . (15) 

As above we take m E to be determined by perturbative matching to order g 2 and g when cal- 
culating the pressure to two and three loop order, respectively, and we neglect Xe because, as we 
have seen above, it starts to contribute only at order g® with small numerical effects. 



A. Two-loop order 

In DRSPT to two- loop order, the result for the pressure is given by the sum of @ and (JSJ) 
except that the latter now involves 

/« +(2) = -f m 3 + A ( I + 41n + 3 ) g 2 E m 2 + ™8m 2 + 5/£ RSPT (16) 
Sir 8tt z \e 2m J ir 

where m 2 = m 2 E + 5m 2 according to ()15|) , The term ^5m 2 arises from a one-loop diagram where 
the counterterm 5m 2 , which itself has to be counted as a one-loop quantity, has been inserted. 

SPT generally produces additional UV divergences and associated scheme dependences, which 
can be seen here in that the pole term in the second term on the right-hand side of Eq. Q16|) no 
longer matches (J7J). Following Refs. 11, 2> Q E3- E3 we treat those by minimal subtraction. 
This means that we introduce a counterterm <5/^ RSPT with 5/^ RSPT — 5/e oc g E 5m 2 j, which can 
be discarded in the perturbative matching as being of higher order and will disappear in fact at 
the next loop order (see below). However, the replacement of m 2 E by m 2 E + 5m 2 has the effect of 
modifying the dependence on at the present loop order. 

The untruncated two-loop result, including now ()16|) in place of @, can be optimized by a 
variational principle (principle of minimal sensitivity) for 5m 2 . This leads to the variational mass 



m 2 = m\ -3a s (p)Trn ( 41n ^ + 1 ) • ( 17 > 



As before and similar to what is done in Refs. |17l.ll8l] we equate A^; and fj,. Since the latter is always 



set proportional to T, this give s rise to a term 5m 2 oc g 3 T 2 ln(g), like in the actual next-to- leading 
order Debye mass in QCD [38j |. however both coefficient and sign are different (and the logarithm 
involves an UV rather than an IR cutoff). This is no contradiction, however, since the variational 
mass (like the effective mass parameter m£) is not identical with the Debye mass responsible for 
exponential screening behavior of the full electrostatic propagator. In fact, we shall see below that 
at higher loop orders the gap equation does not contain a logarithmic term. 
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The simple "gap equation" (|17|) is different from the one obtained in 2-loop HTLPT [17 . 
It is much simpler in form, and it turns out to have a numerical solution that behaves differently 
from that of 2-loop HTLPT. While the latter increases sharply at large coupling, the solution of 
(|17|) saturates at the value m ~ 0.503T as a s increases when Nf = and A# = p, = 2irT [for 
larger /smaller fl the saturation occurs at smaller /larger values; for ft < lAirT there is a maximum 



11 When fermions are included in 2-loop HTLPT [T^| the latter gives rise to a gap equation for fermions as well, 
whereas in DRSPT fermions contribute only through the parameter tue- 
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FIG. 6: The two-loop DRSPT result for pure-glue QCD (gray band) in comparison with the two-loop 
HTLPT result of Ref. 17] (dashed lines) and the "NLA" result of Ref. 20] (dash-dotted lines), all with p, 
varied around 2-kT by a factor of 2. 

value of a s beyond which solutions no longer connect continuously to the perturbative leading- 
order result — for instance, for our lowest value p = ttT this restricts a s to less than 0.34, which 
however presents no problem to the following application]. 

The thermal pressure at two-loop DRSPT is obtained by evaluating (|16j) at the solution of (|17j) 
and combining it with the hard contribution ijljl. The result for pure-glue QCD is given in Fig. |2 
where it is compared with the untruncated two-loop result which uses the perturbative leading- 
order mass instead of the solution of the gap equation (|17f) . In Fig.|f)]the two-loop DRSPT result is 
furthermore c omp ared to the result of the two-loop HTLPT calculation of Ref. jrjj and the "NLA" 
result of Ref. [2TJ which is based on the ^-derivable two-loop expression for the entropy evaluated 
with HTL propagators that include next-to- leading order corrections to asymptotic thermal masses. 
We find that the two-loop DRSPT result has a rather small scale dependence like the two-loop 
HTLPT result, but is significantly below the latter and thus closer to the lattice data. The nonlinear 
scale dependence of the DRSPT result in fact makes it possible to eliminate the scale dependence 
by a principle of minimal sensitivity. The result is given by the lower boundary of the DRSPT 
band in Fig. El This optimized DRSPT result happens to lie right at the center of the estimated 
error of the NLA result of Ref. 0] for T > 3T C [the fact that the latter sharply drops close to 
T c is in fact not a prediction of the NLA result but comes from the necessity to fix an integration 
constant which has been chosen such that P(T C ) = 0]. 

Here we have neglected the soft 2-loop contribution involving Xe of Eq. Q because HTLPT 
does not include a comparable term. When (JSJ) is included, the 2-loop gap equation is modified 
by an extra contribution — ^(9 — Nf)a 2 s Tm on the right-hand side of Eq. (|17|) . Its (rather small) 
effect is displayed on the occasion of the comparison with 4-loop DRSPT in Fig. [Tj 



B. Three-loop order 



The three-loop free energy of EQCD in DRSPT can be easily derived from the results of Ref. 
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+ ™ § m 2 + ±- g %5m 2 - -^-(5m 2 ) 2 . (18) 
7T Air 2 ^ Anm 

where m? = m 2 E + 5m 2 is now defined with the order- g 4 result for m E , Eq. (jlUj) . and where up to 
two insertions of the SPT counterterm 5m 2 have to be included. 

At this order the 5m 2 counterterm of DRSPT restores the UV divergent part to be proportional 
to m 2 E (second term on the r.h.s. of Eq. (|T5|l). and because no further pole terms appear at three- 
loop order, the additional UV divergences of SPT drop out altogether. 

The attempt to determine m by a variational principle leads to the simple quartic gap equation 

G®> (m) S \{m 2 - m 2 E ) 2 + ^ (m 2 -m 2 E )m + (^) ' (g - H ln 2 + ^ m2 = , (19) 

which we have written out for general color number N c . However, one can readily prove that this 
equation has no solution that connects continuously to the perturbative result mg as a s — ► (for 
any value of N c ; also inclusion of the term Q does not change this situation). 

The same problem in fact occurs in SPT in scalar cf) 4 theory to three-loop order Like 
Ref. |.ll| one may look for alternative prescriptions to set up a gap equation. The simplest con- 
ceivable option, however, clearly is to set dm 2 = 0. This trivially connects to perturbation theory, 
and amounts to our previous suggestion of keeping the 3-loop contributions of EQCD untruncated, 
which is in fact most natural since rriE in EQCD is just a mass parameter. Fortunately, the now 
nonlinear scale dependences can be eliminated by a principle of minimal sensitivity, as we have 
discussed above, so in a sense our above improvement may be taken as a trivial implementation 
of DRSPT, with the variational parameter being the renormalization and separation scale rather 
than an additional mass term. 



C. Four-loop order 



In DRSPT to four-loop order (neglecting Xe contributions) we have 
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(20) 



which is to be added to (|18jl . evaluated now also with g 2 E to order g 4 as given by (|13|) . 

The variational gap equation associated with the sum of (|18|) and (|2Uj) does not involve the 
unknown constant c and is therefore known completely to order g E (in fact, the contribution (JSJ) 
proportional to does not enter either, while \ 2 E contributions to (|2*U)) are already of order g 9 or 
higher). This equation is now of sixth order in m and reads 



G (4)( 
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(m 2 



m 2 E 



)G (3) (m) + 2 
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(21) 



As it turns out, the inclusion of the 4- loop logarithm lifts the impasse encountered with DRSPT 
at 3-loop order: there is now a (unique) solution to the 4-loop gap equation which does connect 
continuously to perturbation theory and which is given by the quadratic gap equation 12 



m 



m\ 



0.33808 g z E N c m, 



(22) 



12 The numerical coefficient therein is given by the (real) root of a cubic equation involving the somewhat unwieldy 
constants appearing in If 811 and and could in principle be given in closed (but lengthy) form. 
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FIG. 7: The 4-loop DRSPT result for pure-glue QCD (3 full lines corresponding to p, = ttT, 2ttT, 4ttT and 
5 = 1/3) in comparison with the two-loop DRSPT result (now evaluated with the contribution JHJ included) 
and the untruncated 4-loop with S — 1/3. 

which appears as a factor of (m) and whose only solution with real and positive m is given by 

m = yJm 2 E + (0.16904 g 2 E N c ) 2 - 0.16904 g%N c . (23) 

It is intriguing that the quadratic gap equation (I22|) is of the same form as the one adopted in 
the NLA approximation of Ref. 0] for the asymptotic thermal masses, and also the coefficients 
in the two gap equations happen to be very close (0.338iV c versus N c /ir). However, it should be 
emphasized that the gap equation of DRSPT has no physical meaning outside of DRSPT. Indeed 
there is no reason to expect the leading correction to tre as prescribed by the 4-loop DRSPT gap 
equation to remain the same at 5-loop level, if at that order solutions exist at all. At any given loop 
order, the deviation of m from rriE only influences the orders beyond the perturbative accuracy. 

Numerically, the deviation of m as given by (|22j) from the perturbative value the has some 
effect as displayed in Fig. [Tj It turns out that the 4-loop DRSPT differs from the untruncated 
4-loop result of the previous section in that it has a significantly larger g 7 coefficient (by a factor 
of almost 6 when Jl = 2ttT). As a consequence, this gives a slightly different "prediction" for the 
unknown g G constant 5, but otherwise the results are quite similar to those obtained in the simple 
untruncated evaluation. 

In fact, at 4-loop order the order-g 7 coefficient could in principle be calculated completely, if m E 
is determined to 3-loop accuracy and relevant higher-dimension operators in the effective theory 
are included. (The order g 6 coefficient of course remains beyond the reach of perturbation theory.) 
In this case, however, 4-loop DRSPT would spoil the achieved perturbative accuracy, because it 
changes the g 7 coefficient without correcting these changes through the SPT counterterms, which 
only take care of orders up to and including m 5 T[(5m 2 ) / m 2 ] 3 ~ g®T 4 . Thus, beginning at 4-loop 
order, (DR)SPT ceases to be a possible improvement over (truncated or untruncated) perturbation 
theory. 
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FIG. 8: The exact result of the thermal pressure in the limit of large Nf from Ref. [4fJ|, normalized to 
that of free gluons and as a function of g 2 Nf(ft — 7re~ 7 T) (full line), in comparison with the untruncated 
two- loop and three- loop results (darker and lighter gray areas, respectively), each with ft varied around a 
central value of 2irT by a factor of e. The dash-dotted line corresponds to fastest apparent convergence. 



IV. LARGE- TV/ LIMIT 



We finally consider also the recently solved large-iVf limit of QCD 3{J 4(3, which has been 
proposed as a testing ground for improvements of perturbation theory. In this limit only terms 
involving products a s Nf are kept in the above results and a s itself is taken to zero. The dimen- 
sionally reduced theory is therefore non-interacting in the large- Nf limit and the soft contributions 
are given exactly by ju = ~^ m %- Still, at a given loop order for the hard contributions, we can 
investigate the difference between a strictly perturbative evaluation of /m versus an untruncated 
one which resums an infinite number of terms with odd powers in g. Also the gap equations of 
DRSPT become trivial (but solvable): they all amount to setting m? = m E . 

The results of a numerical evaluation of the untruncated two-loop and three-loop results are 
compared with the exact result of Ref. j4jj in Fig. |SJ Again one can observe a great reduction 
of the scale dependence by going from two-loop to three-loop order. Compared to the strictly 
perturbative result to order g 5 (not displayed in Fig. |HJ) the reduction of the size of the scheme 
dependence is less important than in the pure-glue case, e.g., at g 2 Nf = 10 the reduction is about 
16%. 

Also in contrast to the pure-glue case, the three-loop result has a monotonic scale dependence, 
so the scale cannot be fixed by minimal sensitivity. As remarked in Ref. 40], \x could instead be 
fixed by fastest apparent convergence (FAC). Requiring e.g. that the a 2 term in m 2 E vanishes leads 
to the choice p,/T = ire^^ 1 « 0.93irT. For this value the untruncated 3-loop result coincides with 
the result to order g 5 in strict perturbation theory, which in turn agrees quite well with the exact 
result up to rather large coupling [40] . So while this comparison does not favor one over the other, 
it shows that with an optimal choice of the renormalization scale both the perturbative result to 
order g 5 and the untruncated 3-loop result fare rather well in the large- Nf limit. 

If one applies the same prescription to the 3-loop result of EQCD in pure-glue QCD, one is lead 
to setting p,/T = 47re _ ' r ~ 5 / 22 ~ 1.79irT. This is lower but close to the scale selected by minimal 
sensitivity of the untruncated 3-loop result. Correspondingly, the numerical result following from 
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the FAC choice is fairly close to the one obtained from minimal sensitivity (dash-dotted and full 
line in Fig. El respectively), at least for T > 3T C , which appears to validate the PMS results. 13 

V. DISCUSSION AND CONCLUSION 

We have found that the scale dependence and convergence of the results for the thermal pressure 
from perturbative QCD at high temperature are significantly improved when the contributions from 
soft scales as given by the effective dimensionally reduced theory EQCD are not treated in strict 
perturbation theory. 

In particular, we have explored the predictions of a simple loop expansion of EQCD, in which, 
after strictly perturbative matching of the parameters of the effective Lagrangian, the results are 
not subsequently expanded out in powers of g and truncated. The result obtained at the two-loop 
level in this scheme includes contributions to order g 4 lng completely while being incomplete to 
order g A , and is such that the pressure no longer exceeds the ideal- gas limit. The scale dependence 
is large, however the three-loop result is within the estimated boundaries of the two-loop result. 
This three-loop result has a smaller scale dependence than that of strict perturbation theory to 
order g 5 , and moreover the scale dependence is nonmonotonic so that it can be eliminated by a 
principle of minimal sensitivity (PMS). The correspondingly optimized result is rather close to 
the lattice data on the continuum limit of pure-glue QCD for T > 3T r . When including 4-loop 
effects, in particular the recently determined <? 6 In g contribution of Ref. [23J, we find that this trend 
continues and, in line with [33U37l|. that it is quite possible that all higher-order contributions add 
up to a very small correction above ~ 3T C . 

We have also considered variationally improved "screened" perturbation theory in the dimen- 
sionally reduced theory (DRSPT), where it is a much simpler, gauge invariant alternative to 
HTLPT (though not extensible to dynamic quantities, as HTLPT in principle is). The result 
for the pure-glue pressure when improved through 2-loop DRSPT turns out to be significantly 
lower than that of HTLPT and for T > 3T C fairly close to the lattice results as well as to the 
results of Ref. 0. 

An obvious advantage of DRSPT over HTLPT is that it does not modify the theory at hard 
momentum scales, where the HTL approximation in general breaks down. The latter continues to 
be a good approximation at hard momentum scales only at soft virtuality. In the entropy-based 
HTL-resummations of Ref. |2o| it turns out that the contributions are predominantly coming from 
hard momenta close to the light-cone. In the HTLPT approach, on the other hand, spurious 
contributions at hard momenta occuring at a given loop order are corrected for by the specific 
counterterms of HTLPT at next loop order, so this may present a problem at low loop orders. 
Our two-loop DRSPT results seem to indicate that this is indeed the reason for the difficulties of 
two-loop HTLPT 0, EF 



An unsatisfactory feature of DRSPT, as observed before in SPT applied to scalar field theories 
and thus not unlikely to affect HTLPT as well, is that at three-loop order the mass gap equation 
does not have solutions which connect to perturbation theory. This impasse happens to be lifted 
by the inclusion of soft four-loop logarithms, and the result is then close to that obtained by a 
simple untruncated evaluation of all soft contribtions. Nevertheless, the gap equations of (DR)SPT 



13 The upper boundary of the range of the three- loop results shown in Fig. |H1 follows the exact result to even much 
larger values. There is in fact a choice of p, « 0.75ttT for both the untruncated three-loop result and the strictly 
perturbative result, where they become almost indistinguishable from the exact result. A renormalization scale 
close to this turns out to be favoured by applying FAC to the perturbative result to order g 6 , which has been 
extracted numerically in the large- TV/ limit |4l|. but which is strictly non-perturbative in real QCD. 
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have no particular physical interpretation (as discussed after Eqs. (|17|) and (|23|)). which casts some 
doubt on the systematics of SPT and its usefulness in improving perturbation theory. 

Evidently, our main result is that the convergence behavior of successive approximations to 
the pressure is dramatically improved by abandoning strict perturbation theory in the soft sector. 
Treating this sector beyond strict perturbation theory is in fact closer in spirit to the so-called 
^-derivable app roximations 42] which are the basis for the resummation techniques developed 
in Ref. [l^ . l2(l| |. Such approximations, when implemented in the soft sector, may represent an 
interesting alternative to DRSPT. DRSPT has a single variational parameter, the mass of elec- 
trostatic gluons. While this has the advantage of great simplicity as well as gauge invariance, the 
full self-energy of electrostatic gluons is a nonlocal quantity. One might consider a ^-derivable 
approach which does not have the need for the specific counterterms of SPT and try to construct 
improved approximations, which are in principle gauge dependent, though such gauge dependences 
are strongly suppressed at the variational point |43|. 

It would be interesting to compare such a dimensionally reduced ^-derivable approach with the 
one based on the entropy formalism of Ref. [lii . Eol ]. In the latter, the emphasis is on dynamical 
quasiparticles, which at two-loop order are interaction free. It should be noted that this approach 
is dependent on a real-time formalism which does not lend itself to dimensional reduction. In- 
deed, it involves differentiating thermal distribution functions at the stationary point, where the 
temperature-dependence of spectral functions drops out. However, the relevant theory for the soft 
modes (including hard ones with soft virtuality) is known: at leading order this involves the non- 
local hard-thermal-loop effective action S EU) In fact, the real-time approach might have 
advantages when it comes to including the effects of high chemical potentials /i. In this case dimen- 
sional reduction does not occur. The quasiparticle approach on the other hand appears promising 
for covering the thermodynamics in the entire T-//-plane 42, 48l. I^. 

We intend to investigate these matters in future work. From the present study we conclude 
that perturbative QCD at high temperature is not limited to T > 10 5 T C as previously thought 
0, |(| but, when supplemented by appropriate resummation techniques for soft physics, seems to 
be capable of remarkably good quantitative predictions down to T ~ 2.5T C . 
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